Predict-prevent control method for perturbed excitable systems 
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We present a control method based on two steps: prediction and prevention. For prediction we 
use the anticipated synchronization scheme, considering unidirectional coupling between excitable 
systems in a master-slave configuration. The master is the perturbed system to be controlled, 
meanwhile the slave is an auxiliary system which is used to predict the master's behavior. We 
demonstrate theoretically and experimentally that an efficient control may be achieved. 
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A classical problem in engineering science is the control 
of dynamical systems While earliest applications in 
this field aimed to inhibit the instabilities in electronic 
devices, control studies turned, later on, into other areas. 
Of particular interest are the applications to medicine 
since transitions from regular to irregular oscillations of 
many organs of the human body have been found to be 
associated with a diseased behavior. As an example, the 
heart beat has been found to undergo chaotic behavior 
during arrhythmias and other heart diseases 0. 

One of the best known methods for stabilizing unsta- 
ble periodic orbits is the one introduced by Ott, Grebogi 
and Yorke (OGY) It uses a small perturbation in 
some parameter of the system to stabilize an unstable 
orbit into a periodic one. The delayed-feedback control 
(DFC), introduced by Pyragas j4j, intends the same but 
uses a self-feedback loop and works when the feedback 
delay is close to the period of the unstable orbit that 
one aims to stabilize. It has been experimentally im- 
plemented in optical systems chaotic flows @, and 
cardiac systems 0], among others. 

In this paper we present a rather different control 
method with particular application, but not necessarily 
limited, to excitable systems. These systems have been 
used to model the behavior of many cell types, including 
the typical spiking dynamics of heart cells and neurons. 
Excitable systems are characterized by a highly non- 
linear response to an external perturbation: if the pertur- 
bation is sufficiently small, the system returns smoothly 
to its steady state; if the perturbation exceeds a certain 
threshold, the system fires an excitable spike, sometimes 
called the action potential. During a refractory period 
following the spike, perturbations of moderate amplitude 
do not alter much the dynamics of the system. Although 
excitable systems can involve a large number of variables 
(as in the case of neurons), their essential features can 
be captured with a much-reduced description. We will 
be using in this paper the one-dimensional Adler and 
the two- variable FitzHugh-Nagumo systems as examples. 
They both provide some of the simplest representation of 



excitable dynamics. 

The proposed control method relies on the phe- 
nomenon of anticipated synchronization. Synchroniza- 
tion refers to the collective timing of coupled systems, 
from simple oscillators and clocks to some vital functions 
found in living organisms, among many others [H,0]. An- 
ticipated synchronization refers to a particular regime, 
first studied by Voss [1(3], which appears in unidirection- 
ally coupled systems in a master-slave configuration. In 
this regime, two dynamical systems synchronize in such 
a way that the slave system, y(t), anticipates the trajec- 
tory of the master, x(t). One of the simplest schemes is 
based upon the general equations: 



x(t) = f[x(t)], 

y(t) = f[y(*)] 



K[x(t)-y(t-r)]. 



(1) 
(2) 



The function f[x] defines the dynamical system under 
consideration; K is a general coupling function satisfy- 
ing K[0] = and r is the delay time in the feedback 
loop of the slave. As proven by Voss [l(J, the manifold 
y(i) =x(i + r), i.e. the slave anticipates by a time r 
the dynamics of the master, can be a (structurally) sta- 
ble solution of equations (JT][3J) . Anticipated synchroniza- 
tion has been studied theoretically and experimentally 
in many systems 
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(see [2l| for a short review). 

Taking advantage of anticipated synchronization we 
propose a control method, which we apply to exter- 
nally perturbed excitable systems, designed to correct 
unwanted dynamical behaviors. The method displays 
some of the most advantageous characteristics of OGY 
and DFC, but nevertheless differs essentially from them. 
The method uses the delayed feedback, as in the DFC 
method, but in an auxiliary slave system to process the 
information and decide whether to activate the control. 
After the decision of activating the control has been 
taken, a small perturbation is applied to the master in 
order to prevent its firing. This procedure resembles the 
perturbative action of OGY but it does not require any 
previous calculations. In the following we explain the 
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will take the values /i = 0.95, K = 0.6, t = 1.5, which 
satisfy the conditions for anticipated synchronization (lBI . 



FIG. 1: Schematic presentation of the predict-prevent control 
method. The input signal I(t) is applied only to the mas- 
ter system x while its influence on the slave y is indirectly 
through the unidirectional coupling Kx. The slave dynamics 
has a delayed loop Ky T = Ky(t — r). The control signal is 
applied whenever a condition for the slave variable is satisfied. 



FIG. 2: The solid curve plots cos[:r(t)], the evolution of the 
master system, while the dotted line is cos[y(t)], the slave sys- 
tem coming from a numerical integration of Eqs. (|3l4[ ). Exter- 
nal perturbations of intensity Jo = 0.64 are applied at times 
t = 200 and t = 250 (small vertical lines at the bottom). 
In the first perturbation (t — 200) we see clearly that the 
slave anticipates the dynamics of the master. In the second 
perturbation (t = 250) the spike of the master has been su- 
pressed by applying a correcting pulse of amplitude e = 0.04 
when the trajectory of the slave crosses the reference level 
cos[xo] = 0.707 (horizontal line). We have used the parame- 
ters fi = 0.95, K = 0.6, r = 1.5 and a reaction time ta — 1. 



new method in some detail. We show how it can be im- 
plemented to suppress unwanted spikes produced by the 
effect of external random perturbations. We prove that 
the control method can work under non-ideal conditions, 
such as the presence of noise in the dynamics and pa- 
rameter mismatch between master and slave. Finally, 
we demonstrate experimentally the efficiency of the pro- 
posed method. 

Although our method is quite general, we consider, 
for the sake of concreteness, two unidirectionally coupled 
Adler systems [22[ in the master-slave configuration: 

±{t) = fi-cos[x(t)]+I(t), (3) 
y(t) = fi-cos[y(t)]+Ksm[x{t)-y{t-T)}. (4) 

I(t) represents an external perturbation and [i, K and 
t are constant parameters. Throughout the paper, we 



2l|. Note that the perturbation is applied only to the 
master system x(t) while it acts upon the slave system 
y{t) only indirectly through the coupling term. We use 
a 2ir periodic function for coupling, appropriate for the 
angular- type variables considered. The slave is driven 
by the master and subject to its own feedback loop, of 
strength K, so producing the anticipated synchronization 
regime. The only modification with respect to Eqs. HJH) 
is the inclusion of the external perturbation I(t) acting 
upon the dynamics of the master. As discussed in [l7j |. 
the feedback term lowers the excitability threshold of the 
slave and allows it to react to the perturbation faster 
than the master. Hence, the slave is used to anticipate 
the future response of the master to this perturbation. 
When an unwanted spike appears in the dynamics of the 
slave before it does in the master, a control mechanism 
is triggered and sent to the master in order to prevent 
that event. The whole scheme is shown in Fig. [T] We 
now give some specific details of the method. 

It is easy to see that the master system x{t) is excitable 
for |yit| < 1 with fixed points at x± = ±arccos(^t), and 
oscillatory for \fi\ > 1. If the system is at the stable 
point x(t) = X- and a perturbation acts upon it for a 
short time, the way in which variable x(t) returns to the 
stable value depends on the magnitude of the perturba- 
tion. Let us assume that we apply a single, delta-like, 
perturbation of magnitude Io at time to > 0, i.e. we take 
I(t) = I S(t — t ), when the system is at the stable fixed 
point x(to) = x_. The effect of such delta function per- 
turbation is to instantly change at to the value of x(t) 
such that x(to ) = x(t^ ) + I . The subsequent dynamics 
depends on the magnitude of Iq ■ If Io < 2 arccos(/Lt) the 
system returns smoothly to the fixed point X- , whereas if 
Io > 2 arccos(^) the variable x(t) executes a full rotation 
(a spike) and returns to the equivalent point 2ir + X-. It 
is the aim of our control method to prevent those spikes. 
We do this by adding to the dynamics of the master, 
Eq.®, a small corrective signal — e5(t — T ) of magnitude 
e at time To- The time To is determined by two factors: 
a correction criterion and a response time. The correc- 
tion criterion determines whether the corrective signal is 
sent, and the response time tji is the time it takes the 
correction action to act after it has been required. 

Let us first consider a simple control method in which 
the correction criterion is that the variable x(t) crosses 
a threshold value xq. This resembles the DFC method 
except that we do not have to fix the delay time to sta- 
bilize a given orbit but the delay will correspond to the 
time it takes to inject a corrective signal into the master 
system when it crosses xq, i.e. whenever x(t) > xo- If 
this happens at time t\, then at time To = t\ + tfi the 
correction acts such that x^q ) = x(T ~) — e, and it will 
be considered effective if x(t) is brought back to the sta- 
ble region. The sooner the correction action is taken, the 
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FIG. 3: Correction of a series of spikes. The external pertur- 
bation I(t) consists of a series of random pulses (vertical lines 
al the bottom of both panels) following a Poisson distribution 
with (t) = 500 (there are 69 excitations altogether). Noise of 
intensity D — 0.01 has been added both to the master and to 
the slave. Top panel: control of magnitude e = 0.3 using the 
slave variable y. Bottom panel: control of magnitude e = 0.85 
using the master variable x. The recovery time is t rec = 5 and 
the response time in = 1, other parameters as in figure [2] 



smaller the magnitude of the correction e needed. In any 
event, even for a zero response time, t R — 0, the condi- 
tion e > xq — x+ must be satisifed in order to bring back 
x{t) from the threshold value xq to the stable region. 

Our method works similarly but the decision of ac- 
tivating the corrective signal depends on the dynamics 
of the auxiliary system y(t), coupled to x(t) using the 
anticipated synchronization scheme discussed previously. 
The main idea is that the dynamics of y(t) anticipates 
that of x(t) and hence the correction criterion y(t) > xa 
happens earlier. Consequently, and this is an important 
advantage of the method, a much smaller magnitude e 
is needed to suppress the undesired firing. Alternatively, 
a larger reaction time can be used to decide whether to 
apply the correcting signal. The correction criterion is 
now modified to the following: if y(t) crosses the thresh- 
old value xq then a delta-like correction of magnitude — e 
is applied to the variable x(t). The same conditions for 
reaction and recovery times discussed before also apply 
here. 

We now compare in detail the control methods using 
the criterions based upon the y(t) or the x(t) variables. 
When using, for example, a threshold value xq = 7r/4, our 
numerical results show that for a response time 4r = 0, 



FIG. 4: Trajectories of the master (upper trace) and slave 
systems (lower trace) in a experimental implementation of the 
anticipated synchronization scheme in an electronic circuit 
simulating the FithHugh-Nagumo neuronal dynamics. 



a value e = 0.468 is needed when using the criterion on 
the x(t) variable, while a much smaller value e = 0.022 
suffices in our scheme using the criterion based on y(t). 
The same result applies in the case of a non-zero response 
time. If tu = 1, e reduces from 0.834, when using the 
criterion of x(t), to 0.0307 when using our method, or 
from 1.751 to 0.0426 if t R = 2. For t R = 3 the method 
based on x(t) is unable to control the system while for our 
method it suffices to take e = 0.0597. In Fig. [2] we show 
that our control methods is efficient and indeed enables 
the suppression of a single spike. 

In a typical control situation, the arrival of the pertur- 
bation and the consequent generation of spikes would be 
random. To model this, we consider a train of pertur- 
bations arriving at times ti such that the time intervals 
U+i ~ U are distributed according to an exponential dis- 
tribution of mean value (t) . We also add to both equa- 
tions §5§ and ((3]) independent white noise terms D£(t) 
with correlations (£(£)£(£')) = S(t — t'). We also consider 
that the control mechanism has a recovery time t rec such 
that t rec must elapse after a control signal has been ap- 
plied before another control signal can be activated. The 
top panel of Fig.[3]shows the result of our control method 
while the bottom panel shows the control method based 
on the x(t) variable. It can be observed that not only 
the intensity e of the control is smaller in the first case, 
but also the spikes are much better suppressed. Finally, 
we have checked that our control method still works well 
when the auxiliary system is not an identical copy of the 
master system. 

In order to assess the practical usage of our method 
we have implemented experimentally an electronic ver- 
sion of the FitzHugh-Nagumo neurons, another type of 
excitable system [23| . The details concerning the experi- 
mental setup can be found in [14[ . For this work the set- 
up was modified introducing a monostable circuit that 
generates a pulse with programable variable width and 
delay when the slave system signal crosses a threshold. 
The pulse is applied as a corrective signal to prevent fir- 
ings of the master. The complete diagram of the set-up 
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imentally the efficiency of the method. The proposed 
technique has been applied to excitable systems but we 
hope that this work will stimulate additional studies in 
order to extend this method to other types of systems, 
specially chaotic ones, as well as to identify possible sit- 
uations of practical interest. 

We acknowledge financial support by the MEC (Spain) 
and FEDER (EU) through projects FIS2006-09966 and 
the EU NoE BioSim LSHB-CT-2004-005137. 
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FIG. 5: Experimental verification of the predict-control 
method in the electronic FithHugh-Nagumo system. The 
small vertical lines at the top of each figure are the external 
perturbations, and the lower traces are the master's variable. 
In the bottom panel, the spikes have been controled using a 
criterion based upon the dynamics of the slave (not shown). 
The top panel shows that the same procedure fails to pre- 
vent the spikes if the criterion is applied using the master's 
variable. 



can be found in 24J. Our experimental findings qualita- 
tively agree with the numerical predictions. Fig. [¥] shows 
the dynamics of the master and the slave in response to 
a perturbation applied only to the first. It can be clearly 
seen that the slave anticipates the firing of the master by 
roughly 400 fis. In Fig.[5]we show time traces taken from 
the oscilloscope. The top panel shows the results when 
the master system is controlled using its own output. In 
this situation it is almost impossible to cancel the firing. 
The bottom panel shows the results using our control 
method based on the anticipating slave. Note that all 
spikes have been suppressed in this particular realization 
of the experiment, a typical situation. 

In conclusion, we have proposed a predict-prevent con- 
trol method for the suppression of spikes in perturbed 
excitable systems. The method consists of two steps: (1) 
prediction made by the slave system with the use of the 
anticipated synchronization scheme and (2) correction 
applied to the master in order to suppress its unwanted 
response. We have demonstrated numerically and exper- 
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